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. Nambu and Jona-Lasinio (Phys. Rev. 124 (1961) 246) is formulated 

from the first principles of QCD. The deuteron appears as a neutron- 
proton collective excitation, i.e. a Cooper np-pair, induced by a phe- 
nomenological local four-nucleon interaction in the nuclear phase of 
QCD. The model describes the deuteron coupled to itself, nucleons 
and other particles through one-nucleon loop exchanges providing a 
"q ■ minimal transfer of nucleon flavours from initial to final nuclear states 

\ and accounting for contributions of nucleon-loop anomalies which are 

completely determined by one-nucleon loop diagrams. The dominance 
of contributions of nucleon-loop anomalies to effective Lagrangians of 
^ , low-energy nuclear interactions is justified in the large Nc expansion, 
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1 Introduction 



In the beginning of sixties Nambu and Jona-Lasinio suggested a dynam- 
ical theory of elementary particles [1,2], the Nambu- Jona-Lasinio (NJL) 
model, in which nucleons and mesons are derived in a unified way from a 
fundamental spinor field on the basis of the relativistic extension of the BCS 
(Bardeen-Cooper-Schrieffer) theory of superconductivity [3]. Nowadays the 
NJL model has found a great support in the form of the extended Nambu- 
Jona-Lasinio (ENJL) model with chiral £7(3) x £7(3) symmetry as an effective 
phenomenological approximation of low-energy Quantum Chromodynamics 
(QCD) [4-6]. Chiral perturbation theory within the ENJL model with a 
linear and a non-linear realization of chiral £7(3) x £7(3) symmetry has been 
developed in Refs. [7,8] and [9], respectively. In the ENJL model mesons are 
described as qq collective excitations, the qq Cooper pairs, induced by phe- 
nomenological local four-quark interactions. In turn, the low-lying octet and 
decuplet of baryons can be considered in the NJL approach as three-quark 
qqq collective excitations produced by phenomenological local six-quark in- 
teractions [8]. As has been shown in Refs. [4-11] the ENJL model with chiral 
£7(3) x £7(3) symmetry describes at the quark level perfectly well strong low- 
energy interactions of hadrons in the form of Effective Chiral Lagrangians 
[12,13]. 

In parallel to the description of mesons as a collective excitations of a 
unified spinor field Nambu and Jona-Lasinio suggested to treat the deuteron 
as a neutron-proton collective excitation, i.e. some kind of a Cooper np-pair 
[2]. A phenomenological local four-nucleon interaction has been written in 
the form [2]: 

Ant(x) = - g [^{x) lfl i) c {x)i) c {x)^i){x) + tP(x)a^ij c (x)ij c (x)a^ij(x) 

+^(x)7 M 7 5 f^ c (x) ■ ip c (x)^ 5 Tip(x)}, (1.1) 

where ip(x) is a doublet of the nucleon field, then ip c (x) = Cip T (x) and 
ijj c (x) = ip T (x) C, C is a charge conjugation matrix and T is a transposition; 
t = (t 1 , r 2 , r 3 ) are the isotopical Pauli matrices, and a^ u = (7^7^ — ■y u j fl )/2. 

According to the Nambu- Jona-Lasinio prescription [2] the phenomeno- 
logical interactions Eq. (|1.1| ) should lead to two bound states: a pseudovector, 
isoscalar (J = 1 + , I = 0) or the deuteron, where J and I are the total spin 
and isospin, respectively, coming from the first two interaction terms, and a 
scalar, isovector (J = + ,I = 1) coming from the last term. Unfortunately, 
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Nambu and Jona-Lasinio did not consider the evaluation of the binding en- 
ergy, the magnetic dipole and electric quadrupole moments of the deuteron 
in their approach to the deuteron as the Copper np-pair. Such attempts 
have been undertaken within the model which has been called the Relativis- 
tic field theory model of the deuteron (RFMD) suggested in Refs. [14,15]. 
The RFMD realizes the consideration of the deuteron in the analogous way 
to the Nambu- Jona-Lasinio approach [2] . For the practical evaluation of the 
low-energy parameters characterizing the deuteron there has been suggested 
in the RFMD that the deuteron couples to itself, nucleons and other particles 
through one-nucleon loop exchanges providing a minimal transfer of nucleon 
flavours from initial to final nuclear states and accounting for contributions 
of nucleon-loop anomalies which are completely determined by one-nucleon 
loop diagrams. Then, there has been argued the dominant role of one- 
nucleon loop anomalies for the one-nucleon loop exchanges describing strong 
low-energy nuclear forces. The main problem of the attempts expounded in 
Refs. [14,15] as well as the original idea of the Nambu and Jona-Lasinio to 
describing the deuteron as the Cooper np-pair [2] is in a poor relation to 
QCD. 

In this paper we show that the consideration of the deuteron as a Cooper 
np-pair induced by a phenomenological local four-nucleon interaction and 
the description of low-energy couplings of the deuteron to itself, nucleons 
and other particles through one-nucleon loop exchanges, where nucleon-loop 
anomalies play the dominant role, is fully motivated by low-energy QCD. 
The deuteron appears as a Cooper np-pair in the nuclear phase of QCD and 
couples to itself and other particles through the one-nucleon loop exchanges. 
The dominance of nucleon-loop anomalies occurs naturally as a consequence 
of the large Nc expansion, where Nc is the number of quark colour degrees 
of freedom [16,17]. Nowadays the large N c expansion suggested by 't Hooft 
[16] is accepted as a non-perturbative approach of low-energy QCD to the 
analysis of strong couplings of hadrons and nuclei at low energies [17]. 

Below we would call the Nambu- Jona-Lasinio approach to the deuteron 
based on phenomenological local four-nucleon interactions like those given 
by Eq.(L~T) as the Nambu- Jona-Lasinio model of light nuclei or shortly the 
nuclear Nambu- Jona-Lasinio model with the abbreviation the NNJL model. 

The paper is organized as follows. In Section 2 we discuss the non- 
perturbative phases of QCD and formulate the NNJL model from the first 
principles of QCD. In Section 3 we derive the effective Lagrangian for the 
free deuteron field induced in the nuclear phase of QCD as the neutron- 
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proton collective excitation (the Cooper np-pair) by a phenomenological 
local four-nucleon interaction. We demonstrate the dominant role of one- 
nucleon loop anomalies for the formation of the effective Lagrangian of the 
free deuteron field. In Section 4 we investigate the electromagnetic proper- 
ties of the deuteron and derive the effective Lagrangian of the deuteron field 
coupled to an external electromagnetic field through the magnetic dipole and 
electric quadrupole moments. We show that the effective Lagrangian of the 
electromagnetic interactions of the deuteron calculated in the one-nucleon 
loop approximation at leading order in the large Nc expansion are defined 
by the anomalies of one-nucleon loop diagrams and have the form of well- 
known phenomenological electromagnetic interactions introduced by Corben 
and Schwinger [18] and Aronson [19] for the description of charged vector 
boson fields coupled to an external electromagnetic field. In the Conclusion 
we discuss the obtained results. 

2 Non-perturbative phases of QCD 

The derivation of the NNJL model from the first principles of QCD goes 
through three non-perturbative phases of the quark-gluon system. We call 
them as: 1) the low-energy quark-gluon phase (low-energy QCD), 2) the 
hadronic phase and 3) the nuclear phase. 

The low-energy quark-gluon phase of QCD can be obtained by integrat- 
ing over fluctuations of quark and gluon fields at energies above the scale of 
spontaneous breaking of chiral symmetry (SB%S) A x ~ 1 GeV. This results 
in an effective field theory, low-energy QCD, describing strong low-energy in- 
teractions of quarks and gluons. The low-energy quark-gluon phase of QCD 
characterizes itself by the appearance of low-energy gluon-field configura- 
tions leading to electric colour fluxes responsible for formation of a linearly 
rising interquark potential. The former realizes quark confinement. 

For the transition to the hadronic phase of QCD one should, first, in- 
tegrate out low-energy gluon degrees of freedom. Integrating over gluon 
degrees of freedom fluctuating around low-energy gluon-field configurations 
responsible for formation of a linearly rising interquark potential one arrives 
at an effective field theory containing only quark (q) and anti-quark (q) de- 
grees of freedom. This effective field theory describes strong interactions at 
energies below the SB^S scale A x ~ 1 GeV. The resultant quark system 
possesses both a chirally invariant and chirally broken phase. In the chirally 
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invariant phase the effective Lagrangian of the quark system is invariant un- 
der chiral U(3) x U (3) group. The chirally invariant phase of the quark system 
is unstable and the transition to the chirally broken phase is advantageous. 
The chirally broken phase characterizes itself by three non-perturbative phe- 
nomena: SB^S, hadronization (creation of bound quark states with quantum 
numbers of mesons qq, qqqq, baryons qqq and so on) and confinement. The 
transition to the chirally broken phase caused by SB%S accompanies itself 
with hadronization. Due to quark confinement all observed bound quark 
states should be colourless. As gluon degrees of freedom are integrated out, 
in such an effective field theory the entire variety of strong low-energy inter- 
actions of hadrons at energies below the SB%S scale A x ~ 1 GeV is described 
by quark-loop exchanges. 

Since nowadays in continuum space-time formulation of QCD the inte- 
gration over low-energy gluon-field configurations can be hardly performed 
explicitly, phenomenological approximations of this integration represented 
by different effective quark models with chiral U(3) x U (3) symmetry moti- 
vated by QCD are welcomed. 

The most interesting effective quark model allowing to describe analyti- 
cally both SB^S and bosonization (creation of bound qq states with quantum 
numbers of observed low-lying mesons) is the extended Nambu-Jona-Lasinio 
(ENJL) model [4-11] with linear [7,8,10] and non-linear [9,11] realization of 
chiral U(3) x U(3) symmetry. As has been shown in Ref. [6] the ENJL model 
is fully motivated by low-energy multi-colour QCD with a linearly rising 
interquark potential and Nc quark colour degrees of freedom at Nc — > oo. 
In the ENJL model mesons are described as qq collective excitations (the 
Cooper gg-pairs) induced by phenomenological local four-quark interactions. 
Through one-constituent quark-loop exchanges the Copper gg-pairs acquire 
the properties of the observed low-lying mesons such as 7r(140), K(A98), 
77(550) , p(770), cj(780), if* (890) and so on. For the description of low-lying 
octet and decuplet of baryons the ENJL model has been extended by the 
inclusion of local six-quark interactions responsible for creation of baryons 
as ggg collective excitations [8]. 

Integrating then out low-energy quark-field fluctuations, that can be 
performed in terms of constituent quark-loop exchanges, one arrives at the 
hadronic phase of QCD containing only local meson and baryon fields. The 
couplings of low-lying mesons and baryons are described by Effective Chiral 
Lagrangians with chiral U(3) x U(3) symmetry [4-13]. 

The nuclear phase of QCD characterizes itself by the appearance of bound 
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nucleoli states - nuclei. In order to arrive at the nuclear phase of QCD we sug- 
gest to start with the hadronic phase of QCD and integrate out heavy hadron 
degrees of freedom, i.e. all heavy baryon degrees of freedom with masses heav- 
ier than masses of low-lying octet and decuplet of baryons and heavy meson 
degrees of freedom with masses heavier than the SB^S scale A x ~ 1 GeV. 
At low energies the result of the integration over these heavy hadron de- 
grees of freedom can be represented in the form of phenomenological local 
many-nucleon interactions. Following the scenario of the hadronic phase of 
QCD, where hadrons are produced by phenomenological local many-quark 
interactions as many-quark collective excitations, one can assume that some 
of these many-nucleon interactions are responsible for creation of many- 
nucleon collective excitations. These excitations acquire the properties of 
observed bound nucleon states - nuclei through nucleon-loop and low-lying 
meson exchanges. This results in an effective field theory describing nuclei 
and their low-energy interactions in analogy with Effective Chiral Lagrangian 
approaches [12,13]. Chiral perturbation theory can be naturally incorporated 
into this effective field theory of low-energy interactions of nuclei. 

We would like to emphasize that in this scenario of the quantum field 
theoretic formation of nuclei and their low-energy interactions nuclei are 
considered as elementary particles described by local interpolating fields. In 
parallel to the Nambu-Jona-Lasinio approach to light nuclei [2] the repre- 
sentation of nuclei as elementary particles has been suggested by Sakita and 
Goebal [20] and Kim and Primakoff [21] for the description of electromag- 
netic and weak nuclear processes. We develop the quantum field theoretic 
approach to the interpretation of nuclei as elementary particles represented 
by local interpolating fields by starting with QCD. 

In this scenario the deuteron, being the lightest bound nucleon state, 
appears in the nuclear phase of QCD as the neutron-proton collective ex- 
citation (the Cooper np-pair) induced by a phenomenological local four- 
nucleon interaction caused by the contributions of heavy hadron exchanges 
at low energies. Through one-nucleon loop exchanges the Cooper np-pair 
with quantum numbers of the physical deuteron acquires the properties of 
the physical deuteron (i) the binding energy Ed = 2.225 MeV, (ii) the mag- 
netic dipole moment /id = 0.857 /xn, where /in is a nuclear magneton, (iii) 
the electric quadrupole moment Qn = 0.286 fm 2 [22] and so on. 

We would like to emphasize that Sakita and Goebal [20] by treating the 
deuteron as an elementary particle described by a local interpolating field 
Dn(x) have calculated the cross section for the photo-disintegration of the 



6 



deuteron 7 + D — > n + p within the dispersion relation approach. The 
more recent analysis of the same process by using the dispersion relations 
has been carried out by Anisovich and Sadovnikova [23] based on the disper- 
sion relation technique developed by Anisovich et al. [24]. The dispersion 
relation approach as well as the NNJL model is a relativistically covariant 
one. Within the dispersion relation approach one deals with directly the 
amplitudes of the process of the deuteron coupled to other particles keeping 
under the control intermediate states in the form of the pole and branching 
point singularities. The residues at pole singularities are defined by the ef- 
fective coupling constants which as usual in the dispersion relation technique 
are taken from experimental data. Unlike the dispersion relation approach 
in the NNJL model developing the Lagrange approach to nuclear forces we 
focus on the evaluation of the effective coupling constants via the derivation 
of the effective Lagrangians of the deuteron coupled to nucleons and other 
particles at low energies. These effective coupling constants are defined in 
the NNJL model by one-nucleon loop anomalies related to high-energy NN 
fluctuations of virtual nucleon (N) and anti-nucleon (N) fields. Thus, the 
dispersion relation approach to the description of low-energy interactions of 
the deuteron and the NNJL model do not contradict but should complement 
each other. 

In the form of the path integral formulation of QCD the non-perturbative 
phases of QCD can be represented by the following sequence of transforma- 
tions. Let us start with the path integral over the quark q, anti-quark q and 
gluon A fields related to a generating functional of quark and gluon Green 
functions and defined by 

Z = J DqDqDA e l S d ' x ^ & *> A ^ . (2.1) 

Integrating over high-energy quark-gluon fluctuations restricted from below 
by the SB%S scale A x ~ 1 GeV we arrive at the generating functional 

Z = J DqDqDae^ d4x/: ?s B [Q,Q,A + a] (22) 

describing strong low-energy interactions of quarks and gluons in the low- 
energy quark-gluon phase of QCD, where A and a are non-perturbative 
gluon-field configurations responsible for the formation of a linearly rising 
interquark potential providing quark confinement and the gluon-field fluctu- 
ations around these gluon-field configurations. 
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Integrating then over the gluon-field fluctuations a we obtain the gener- 
ating functional 

Z = J DqDqe^ f ^ x ^ eff [^' ^ oca ^ m ulti — q couplings] ^ 3) 

This generating functional describes the effective theory of quarks coupled to 
each other at energies of order A x ~ 1 GeV and less. At the phenomenolog- 
ical level the result of the integration over gluon-field configurations can be 
represented in the form of phenomenological local multi-quark interactions 
responsible for the creation of multi-quark collective excitations. This quark 
system is unstable under SB%S and hadronization. By converting quark de- 
grees of freedom into the hadronic ones or differently hadronizing the quark 
system we arrive at the generating functional given by the path integral over 
hadronic degrees of freedom only 

Z = J DM e DB e DM h DB h e i f & x C ^[Me, B e , M h , B h ] ^ ^ 

where Mg, Bi and M^, Bh are local interpolating fields of mesons and baryons. 
The indices I and h correspond to light hadrons with masses of order of 
1 GeV and less and heavy hadrons with masses much greater than 1 GeV. 
For the practical applications to the description of low-energy couplings of 
light and heavy hadrons the effective Lagrangian C eS [M e , B e , M h , B h ] can be 
approximated by Effective Chiral Lagrangians for light hadrons [12,13] and 
heavy hadrons [25] as well. The generating functional Eq.( |2.4j ) describes 
low-energy interactions of hadrons in the hadronic phase of QCD. 

Integrating over heavy baryon degrees of freedom given by the fields Mh 
and Bh we get the generating functional in the form of the path integral over 
the light hadron degrees of freedom 

Z = J DMiDBi e* 1 f ^ x ^ cS \Mt-> Be, local multi — Bi couplings] ^ 5) 

At low energies the result of the integration over heavy hadron degrees of 
freedom can be represented phenomenologically by local multi-baryon cou- 
plings some of which should be responsible for the creation of multi-baryon 
excitations with quantum numbers of nuclei. In term of the local interpolat- 
ing fields of nuclei the generating functional Eq.(|2.5|) acquires the form 



DM e DB e D D D TI D 3 He D 4 He . 

,Z U r „ r fl/f . R„ Pi 3tt 3tt 4t 



x e 



% J d*x C eS [M E , B e , D, 3 H, 3 He, 4 He, . . .] ? (2 6) 
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where D, 3 H, 3 He and 4 He are the local interpolating fields of the deuteron, 
the triton, the helium-3 and the helium-4, respectively. The ellipses stand 
for a possible contribution of other nuclei. The generating functional Eq. ( |2.6| ) 
describes the nuclear phase of QCD, when nuclei couple to each other and 
light hadrons at low energies. Chiral perturbation theory [7,8] is naturally 
incorporated into this theory. 

3 The deuteron as a Cooper np— pair 

In order to describe the deuteron as a Cooper np-pair we introduce a 
phenomenological local four-nucleon interaction caused by the integration 
over heavy hadron degrees of freedom. First, let us consider the simplest 
form of this local four-nucleon interaction 

Ant(x)=-^4(X)/(X), (3.1) 

where gy is the phenomenological coupling constant of the NNJL model 
[14,15], Mn = 940 Me V is the nucleon mass. We neglect here the electro- 
magnetic mass difference for the neutron and the proton. As has been found 
in [14,15] the coupling constant gy is related to the electric quadrupole mo- 
ment of the deuteron Q B : g\ = 2tt 2 Q d M^ [15]. 

The nucleon current with the quantum numbers of the deuteron is 

defined by [14,15] 

f(x) = -i [p c (x)^n(x) - n c (x)^p(x)}. (3.2) 

Here p(x) and n(x) are the interpolating fields of the proton and the neutron, 
N c (x) = C N T (x) and N c (x) = N T (x)C, where C is a charge conjugation 
matrix and T is a transposition. In terms of the electric quadrupole moment 
of the deuteron the phenomenological local four-nucleon interaction Eq. fl3.1|) 
reads 

C mt (x) = ~n 2 Q D jl(x)f(x). (3.3) 

Now let us discuss the behaviour of the phenomenological coupling constant 
g"y/4M^ from the point of view of the large Nq expansion in QCD with the 
SU(N C ) gauge group at iV<7 — > oo [16,17]. Suppose, for simplicity, that the 
phenomenological four-nucleon interaction Eq.( [3.1|) is caused by exchanges 
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of the scalar / (980) and a (980) mesons being the lightest states among 
heavy hadrons we have integrated out. 

Through a linear realization of chiral U(3) x U(3) symmetry and the 
Goldberger-Treiman relation one can find that the coupling constant of <x- 
mesons g CT NN, the gg-scalar mesons, coupled to the octet of low-lying baryons 
should be of order ^ ct nn ~ 0(yJN c ) at N c — > oo. The scalar mesons / (980) 
and a (980) are most likely four-quark states with qqqq quark structure 
[26,27]. In the limit Nc — > oo such qqqq states are suppressed by a factor 
1/Nc [17]. Thus, an effective coupling constant of low-energy NN interaction 
caused by the qqqq scalar meson exchanges should be of order 0(1/ No) at 
Nc — > oo. By taking into account that in QCD with No — > oo the nucleon 
mass M N is proportional to N c [17], M N = N c M q , where M q ~ 300 MeV is 
the constituent quark mass, we can introduce the nucleon mass Mn in the 
effective coupling constant as a dimensional parameter absorbing the factor 
N c , i.e. g Y /4M^. This is also required by the correct dependence of the 
deuteron mass on No- As a result the phenomeno logical coupling constant 
gy turns out to be of order 0(y/No) at Nc — > oo. 

We should emphasize that one does not need to know too much about 
quark structure of heavy hadron degrees of freedom we have integrated out. 
Without loss of generality one can argue that among the multitude of con- 
tributions caused by the integration over heavy hadron degrees of freedom 
one can always find the required local four-nucleon interaction the effective 
coupling constant of which behaves like 0(1/ No) at Nc — > oo. As we show 
below this behaviour of the coupling constant of the phenomenological four- 
nucleon interaction leads to the deuteron as bound neutron-proton state, and 
it is also consistent with the large Nc dependence of low-energy parameters 
of the physical deuteron [17]. 

The effective Lagrangian of the np-system unstable under creation of the 
Cooper np-pair with quantum numbers of the deuteron is then defined by 

C np (x) = n(x) (iYd„ - M N ) n(x) + p(x) (17^ - M N ) p(x) 



^Jt(x)f(x) + ..., (3.4) 



where ellipses stand for low-energy interactions of the neutron and the proton 
with other fields. 

In order to introduce the interpolating local deuteron field we should 
linearalize the Lagrangian Eq.(|3.4|). Following the procedure described in 
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Refs. [4-11] for the inclusion of local interpolating meson fields in the ENJL 
model we get 

C np (x) -> n(x) (17^ - M N ) n(x)+p(x) (iyfi^ - M N )p(x) 
+M 2 Dl(x)D»(x) + gvjt(x)D*(x) + gyj^D^x) + ..., (3.5) 

where M = 2 Mn and D^(x) is a local interpolating field with quantum 
numbers of the deuteron. 

In order to derive the effective Lagrangian of the physical deuteron field 
we should integrate over nucleon fields in the one-nucleon loop approximation 
[2,14,15]. The one-nucleon loop approximation of low-energy nuclear forces 
allows (i) to transfer nucleon flavours from an initial to a final nuclear state 
by a minimal way and (ii) to take into account contributions of nucleon-loop 
anomalies [28-31], which are fully defined by one-nucleon loop diagrams 
[29-31]. It is well-known that quark-loop anomalies play an important role 
for the correct description of strong low-energy interactions of low-lying 
hadrons [4-13]. We argue the dominant role of nucleon-loop anomalies for 
the correct description of low-energy nuclear forces in nuclear physics. We 
demonstrate below the dominant role of nucleon-loop anomalies by example 
of the evaluation of the effective Lagrangian of the free deuteron field. 

The effective Lagrangian of the free deuteron field evaluated in the one- 
nucleon loop approximation is defined by [14,15]: 

J d 4 xC eS (x) = J d 4 xM%Dl(x)D»(x) 

_J d4x J <fx0, _ lkl . {x _ Xl) D ^ x)Dv{xi) A rp^ i; Q), (3.6) 
where the structure function W- v {k-i\ Q) is given by 

U^Q) = / ^tr{ s 1 s . 7 " 1 tA. (3.7) 

J nh { M N — k — Q — ki M N — k — Q J 

The 4-momentum Q = a k\ is an arbitrary shift of momenta of virtual nu- 
cleon fields with an arbitrary parameter a. According to Refs. [29,31] the 
Q-dependent parts of one-nucleon loop diagrams are related to the anoma- 
lies of these diagrams, thereby, the correct evaluation of the Q-dependence of 
one-nucleon loop diagrams is a great deal of importance in the NNJL model 
stating a dominant role of nucleon-loop anomalies. For the evaluation of the 



11 



Q-dependence of the structure function H^ik^Q) we apply the procedure 
invented by Gertsein and Jackiw [29] (see also [15]): 

} d 

IP"(fci; Q) - n^ffcu 0) = dx — U^iki; xQ) = 

J dx 
o 

d 4 k . d f 1 1 



J dx J —Q*—ti{ — - -V 



o 
i 



A tr<^ x 

dk x \ M N - k - xQ - fci M N -k-xQ 



= 2 J dx Hm /^-tr{(M N + k + xQ + £1)7" (M N + fc + zQh'M = 

o ^ ' 

= 2 (2Q^ - Q 2 gT) + lik^Q? + A:^g M - h ■ Q g^) = 
= -2 a(a + 1) (ifej ^ - 2 Jfe^). (3.8) 

Thus, we obtain 

^"(fcuQ) -^"(fcuO) = -2a(a + l)(A; 2 ^ 1/ -2£;^). (3.9) 

We would like to emphasize that the r.h.s. of Eq.( |3.9| ) is an explicit expression 
completely defined by high-energy (short-distance) NN fluctuations, since 
the virtual momentum k is taken at the limit k — > oo, and related to the 
anomaly of the one-nucleon loop diagram with two vector vertices (the VV- 
diagram) [29,31]. 

The structure function (/ci ; 0) has been evaluated in Refs. [14,15] and 
reads 

W v (h- 0) = -^k\g^ - ^)J 2 (M N ) + 2< ? HJ 1 (M N ) + M 2 J 2 (M N )], (3.10) 

where Ji(Mn) and J 2 (Mn) are the quadratically and logarithmically diver- 
gent integrals [14,15]: 



Ji(M K 



rf 4 fc 1 _ ? d\k\k 2 
^HM^-k 2 ~ J (M 2 + Py/ 



T /,.v /" d 4 fc 1 „ ^ d\k\k 2 /o,,x 

^< A « ^ / ^ = 2 / MTI^' (3 ' n) 

The cut-off Ad restricts from above 3-momenta of low-energy fluctuations 
of virtual neutron and proton fields forming the physical deuteron [14,15]. 
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As has been shown in Refs. [14,15] the cut-off Ad is much iess than the mass 
of the nucleon, Mn Ad [14,15]. This leads to the relation between the 
divergent integrals: 

MM N ) = 2Ml J 2 (M N ) = \ ^ ~ 0(1/N C ) (3.12) 



which we use below. Note that in Eq.( |3.10|) we have taken into account only 



the leading terms in the external momentum expansion, i.e. the /^-expansion 
[14,15]. 

The justification of the dominance of the leading order contributions in 
expansion in powers of external momenta can be provided in the large Nq 
approach to the description of QCD in the non-perturbative regim. Indeed, 
in QCD with the SU(Nc) gauge group at Nc — > oo the baryon mass is 
proportional to the number of quark colours [17]: Mn ~ Nc- Since for 
the derivation of effective Lagrangians describing the deuteron itself and 
amplitudes of processes of low-energy interactions of the deuteron coupled 
to other particles all external momenta of interacting particles should be 
kept off-mass shell, the masses of virtual nucleon fields taken at Nc — > oo 
are larger compared with external momenta. By expanding one-nucleon 
loop diagrams in powers of 1/Mn we get an expansion in powers of 1/Nq- 
Keeping the leading order in the large Nc expansion we are leaving with the 
leading order contributions in an external momentum expansion. We should 
emphasize that anomalous contributions of one-nucleon loop diagrams are 
defined by the least powers of an external momentum expansion. This implies 
that in the NNJL model effective Lagrangians of low-energy interactions are 
completely determined by contributions of one-nucleon loop anomalies. The 
divergent contributions having the same order in momentum expansion are 
negligible small compared with the anomalous ones due to the inequality 
Mn Ad and the limit Nc — > oo. This justifies the application of the 
approximation by the leading powers in an external momentum expansion to 
the evaluation of the effective Lagrangians of the deuteron coupled to itself 
and other fields. 

Collecting all pieces we get the structure function rP"(fci; Q) in the form 

U^(h;Q) = ^(klg^-k^)J 2 (M N )+2g^[J 1 (M N )+M^J 2 (M N )} 

-2a{a + \)(k\g» v -2kVkX). (3.13) 
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The effective Lagrangian of the free deuteron field is then defined by 
C eS (x) 



x) 



+ M 2 



9v_ 

2tt 2 



[Jx(M N ) + Ml J 2 (M N )] 



(3.14) 



where D^ v (x) = d^D v (x) — d u D fl (x). We have dropped some contributions 
proportional to the total divergence of the deuteron field, since d lJL D^{x) = 0. 
For the derivation of Eq. (|3.14|) we have used the relation 



^ J ^^-i kl -(x- Xl ) D ^ x)Dv(xi) ^_ m 
l -J d'xDl v {x)D^{x). 



(3.15) 



In order to get a correct kinetic term of the free deuteron field in the effective 



Lagrangian Eq. (|3.14|) we should set 



2tt 2 



a(a + 1) = 1. 



(3.16) 



Since a is an arbitrary real parameter, the relation Eq. fl3.16| ) is valid in the 
case of the existence of real roots. For the existence of real roots of Eq. fl3.16|) 
the coupling constant gy should obey the constraint g\ < 8ir 2 that is satisfied 
by the numerical value gy = 11.319 calculated at N c = 3 [15]. Since gy ~ 
0(y/Nc) at Nc — > oo, Eq. (|3.16|) has real solutions for any Nc > 3. 

Due to Eq.( p,16| ) the effective Lagrangian of the free deuteron field takes 
the form 



+ Ml 



x 



9v_ 
2vr 2 



[Jx(M N ) + Ml J 2 (M N )] DUx)D^(x). (3.17) 



By performing the renormalization of the wave function of the deuteron field 
[14,15] 



1 + ^J 2 (M N ) 



1/2 



DHx) -> DHx) 



(3.18) 



14 



and taking into account that Mn 3> we arrive at the effective Lagrangian 
of the free physical deuteron field 

£ eS (x) = ~ Dl(x)D^(x) + Ml Dl(x)D»(x), (3.19) 

where Mq = M — e-d is the mass of the physical deuteron field. The binding 
energy of the deuteron e D reads 

17 g Y J t (M N ) 17 3 n(UN , ,o 9m 

£d = 48 ~mr = is 0(1/iVc) - (3 - 20) 



We have used here the relation between divergent integrals Eq. fl3.ll ) and 



expressed the phenomenological coupling constant gy in terms of the electric 
quadrupole moment of the deuteron g Y = 2ii 2 Q D M^. The dependence of the 
physical observable parameter, the binding energy of the deuteron, on the 
cut-off Ad is usual for any effective theory like the NJL model [1-11]. 

At Nc — > oo the binding energy of the deuteron behaves like 0(1/Nc) 
as well as the electric quadrupole moment Qd and the coupling constant of 
the phenomenological local four-nucleon interaction Eq.( [3.1|) . This testifies 
a self-consistency of our approach. Really, all parameters of the physical 
deuteron field are of the same order according to the large Nc expansion. 
This means that the vanishing of the coupling constant of the phenomeno- 
logical four-nucleon interaction Eq. (|3.1|) in the limit Nc — > oo entails the 
vanishing of all low-energy parameters of the physical deuteron. 



4 Electromagnetic properties of the deuteron 

The description of the deuteron as a Cooper np-pair changes the analysis 
of the electromagnetic parameters of the deuteron given in Ref. [15], since we 
do not have more a "bare" deuteron field having the magnetic dipole and 
electric quadrupole moment. Therefore, for the Cooper np-pair both the 
magnetic dipole and electric quadrupole moments have to be induced fully 
by the one-nucleon loop contributions. For the self-consistent description 
of the electromagnetic properties of the deuteron we cannot deal with only 
the nucleon current j^(x) given by Eq.fl3.2p and have to introduce the tensor 
nucleon current [14,15] 

J^( x ) = pc( x )a^n(x) - n c (x)a^p(x), (4.1) 
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where = (7^ - 7^) /2. 

The local four-nucleon interaction producing the deuteron as a Cooper 
np-pair reads now 

C b *(x) = --^Jt(x)J*(x). (4.2) 



The baryon current J^(x) is defined by 

J M (a;) = — i gy [p c (x)^n(x) — n c (x)^p(x)] 



,/T d u [p c (x)a^n(x) - n c (x)a utM p(x)\, (4.3) 



M N 

where is a dimensionless phenomenological coupling constant [15]. The 
contribution of the tensor nucleon current looks like the next-to-leading term 
in the long-wavelength expansion[] of an effective low-energy four-nucleon 
interaction. 

The effective Lagrangian of the np-system unstable under creation of the 
Cooper np-pair with quantum numbers of the deuteron is then defined 

£ n P( x ) = n(x) (zYfy - M N )n(x) + p{x){i'fd ll - M N )p(x) 

-^4(x)nx). (4.4) 

The linearalized version of the effective Lagrangian Eq. fl4.4|) containing the 
interpolating local deuteron field reads 

C np (x) -»■ n(x)(iYd„- M N )n(x) + p(x)(i^d fl - M N )p(x) 

+M 2 Dl(x)D^x) + g Y ji(x)D»(x) + g v f{x)Dl{x) 

+^Jl(x)LT(x) + §- o J^(x)Dl(x), (4.5) 

where Mq = 2Mn, D^(x) is a local interpolating field with quantum numbers 
of the deuteron and D» v {x) = d»D u (x) - d u D»{x). 

The interactions with the tensor current give the contributions only to 
the divergent part of the effective Lagrangian of the free deuteron field de- 
termined now by 

= -~ ( - ^(a + 1) + ^ + 6 ^f 2 T + ^ UM^MD^x) 
+ (m 2 - ^[MM^+M^MM^D^D^x). (4.6) 



1 Duc to proportionality Mn ~ Nc this expansion is related to the large Nc expansion. 
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The one-nucleon loop diagrams defining the effective Lagrangian Eq.(O) 
are depicted in Fig.l. 

Due to the relation Eq. ( p,16| ) the effective Lagrangian of the free deuteron 
field Eq. (|4.6|) takes the form 

£-(*) = »v + 68vST + 3 g ^ 2(MN) j^ M ^ M 

+ ( M ° " ^5 |Ji(Mn) + KMM s )]\Dl(x)D'(x). (4.7) 

After the renormalization of the wave function of the deuteron field we arrive 
at the effective Lagrangian defined by Eq. (|3.19| ) with the binding energy of 
the deuteron depending on g v an d <7t 

17<fl Ji(M N ) / 48 g T 24g*\ 
£D " 48^ M N ^ + 17 9v 17^j" 



18 QdA H 17^ 17^'' ( 8) 



where we have used the relation between divergent integrals Eq. (|3.11| ) and 
expressed the phenomenological coupling constant g\ in terms of the electric 
quadrupole moment of the deuteron g\ = 2tt 2 QdM^. In order to make 
the prediction for the binding energy much more definite we have to know 
the relation between the phenomenological coupling constants gy and g^. 
For this aim we suggest to consider the electromagnetic properties of the 
deuteron. 

Including the electromagnetic field by a minimal way <9 M — > d^ + i e A^x), 
where e and A^x) are the electric charge of the proton and the electromag- 
netic potential we bring up the linearalized version of the Lagrangian Eq. ( |f.5|) 
to the form 

C np (x) — > £ E lm( x ) = 

= n(x)(i'fdp - M N )n(x)+p(x)(iYd ll - M N )p(x) + M^Dl(x)D^(x) 
+g v jl(x)D»(x) + gyf(x)Dl(x) + J^^(x) + ^J^D\ v [x) 

- ep(x)Yp(x) A, (x) - ie ^Jl(x)(A^x)D^x) - A»(x)D»(x)) 

+ie^-J^(x)(A tl (x)D 1/ (x) -A"(x)D»(x)). (4.9) 
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By using this Lagrangian we should calculate fully all contributions to the 
effective Lagrangian of the deuteron coupled to an external electromagnetic 
field. These are the effective Lagrangians of the Corben-Schwinger [18] and 
the Aronson [19] type defining at field theoretic level the magnetic dipole 
and the electric quadrupole moment of the deuteron, and the effective inter- 
actions which can be identified with the contributions caused by the minimal 
inclusion of the electromagnetic field d fl D u (x) — > (d^ + ie A^(x))D v (x). 



4.1 The phenomenological Corben— Schwinger interac- 
tion 

The one-nucleon loop diagrams defining in the NNJL model effective 
electromagnetic interactions of the deuteron linear in electric charge e in- 
duced by the Lagrangian Eq. (|4.9| ) are depicted in Fig. 2. One can show that 
in the 1/Mn expansion corresponding to the large Nc expansion due to the 
proportionality Mn ~ Nc [17] the one-nucleon loop diagrams in Fig. 2a and 
2b are divergent. Therefore, due to Eq. fl3.12p at leading order in the large 
Nc expansion the contributions of these diagrams can be neglected with 
respect to the contributions of the diagrams in Fig. 2c and 2d defining the 
phenomenological Lagrangians of the Corben-Schwinger, £qs{x), and the 
Aronson, Ca{x), type, respectively, in terms of the nucleon-loop anomalies 
[15]. 

The effective Lagrangian of the diagram in Fig. 2c is defined by [15] 

J d A x C Fig . 2c (x) = J d A x J - d ^t)* 2 D ^ X ^ A v( x z) 

x e - ikvxi e ~ ik2 - X2 e i[kl+k2) - x ^ J^(h, k 2 ; Q). (4.10) 

Att 2 

In the one-nucleon loop approximation the structure function J^^iki, k 2 ; Q) 
is given by the momentum integral [15] 

d A k 



J^(h,k 2 ;Q)= J 



ir 2 i 



xtr^ s ^ a s s s s s (4.11) 

1 M^-k-Q Mn-k-Q-h M^-k-Q -k 1 -k 2 ) 

The 4- vector Q = ak\ + bk 2 , where a and b are arbitrary parameters, 
displays the dependence of the k integral in ( [4.1 1|) on a shift of a virtual mo- 
mentum. According to Refs. [29,31] a Q-dependent part of an one-nucleon 
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loop diagram is related to the anomaly of this diagram. Therefore, the eval- 
uation of the Q-dependence of the one-nucleon loop diagram should play an 
important role in the NNJL model. For the evaluation of the Q-dependence 
of the structure function Eq.( |4.11|) we apply the method invented by Gertsein 
and Jackiw [29] and consider the following difference of momentum integrals 

5J^(kx,k 2 ;Q) = J^(h,k 2 ;Q) ~ J Pa >*{kiM^) (4-12) 

In accordance with the Gertsein-Jackiw method the difference ( |4 . 1 2| ) can be 
represented by the integral 



8J^(k u k 2 -Q) = f 1 dx^-J^(k u k 2 ;xQ) 
Jo ax 



dx / — — Q ——r tr< 7^ - t - f s 



IT 2 I 
1 



dk x { M N — k — xQ M^-k-xQ-h 



-f r ■ (4.13) 

Mn — k — xQ — ki — k 2 ) 

This shows that the contribution of the Q-dependent part of the structure 
function Eq.( [4.11| ) is just the surface term. Following Gertsein and Jackiw 
[29] and evaluating the integral over k symmetrically we obtain 

5J (3atl (k 1 ,k 2 ;Q) = -2 C dx lim ( ^tr{ 7 ^ (M N + k + xQ) 7 a x 

JO k^oo \ fc 4 

x (M N + k + xQ + h) Y (M N + k + xQ + h + k)}j- (4.14) 

The brackets < . . . > mean the averaging over k directions. Due to the limit 
k — > oo we can neglect all momenta with respect to k. 

5J l3a>l (k 1 ,k 2 ;Q) = -2 lim / ^tr{/ j^fry' 4 k}\ (4.15) 

i?^oo \ fc 4 / 

Averaging over /c-directions 

lim — — = - {g^fo + g^g™ + g^gT) (4.16) 

fc^oo fc 4 24 



we obtain 



SJ^(h,k 2 ;Q) = -^tr(7 A7 V7 a Q7 /i + 7 /3 7A7V07' 1 

+ T/tfVVyV) = I (Q a g^ + Q p g» a + Q"g aP )- (4-17) 
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Our result Eq.( 4.17| ) agrees with the statement by Gertsein and Jackiw [29] 
that the Q-dependence of one-nucleon loop diagrams, i.e. the anomaly of 
the one-nucleon loop diagram, is fully defined by the surface behavior of the 
integrand of the momentum integral at a virtual momentum going to infinity, 
k — > oo. This relates the anomalies of the one-nucleon loop diagrams with 
contributions of high-energy (short-distance) fluctuations of virtual nucleon 
and anti-nucleon fields, i.e. the NN fluctuations. 

Now we can proceed to the evaluation of J^ik^ k 2 ; Q). In order to pick 
up the contribution of the Q-dependent part one cannot apply the Feynman 
method of the evaluation of momentum integrals like ( |4.11| ). This method 
involves the mergence of the factors in the denominator with the subsequent 
shift of a virtual momentum. On this way one can lose the Q-dependence by 
virtue of the shift at the intermediate stage. Thereby, we have to evaluate 
the integral over k without any intermediate shifts. 

One can make this by applying a momentum expansion related to the 
1/Mn expansion or that is the same the large Nc expansion [17] and keeping 
only the leading terms. 



J^{ki,k 2 ;Q) = 



d^k \ q M N + k + Q 
tr<7 P 



TC 2 i 



M? 



N 



X 



X 



1 + 



2k-{Q + h) 



Ml 



k 2 



r 



k 2 



1 + 



2k-Q 



M N + k + Q + h 



M 2 -k 2 _ 



M? 



N 



k 2 



X 



M 2 - k 2 



1 + 
rd A k 



2k-(Q + k 1 + k 2 ) 



M 2 -k 2 



rtr 



{M^(k + Qh a Y + M^ a (k + Q + h)^ 



J 7T 2 l {M 2 -k 2 f 

+M^^ a ^(k + Q + h + k)+ l P {k + Qh a (k + Q + h) 

r x | 2 k ■ (3Q + 2ki + k 2 ) 



X7 



'(k + Q + h + k)} 



M 2 



N 



- k< 



1 rd 4 k 
2J 



Ml 



{Ml - k 2 ) 2 (Ml - k 2 f_ 
+^ a ^(Q + h + k 2 )} 
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+^k 1 a h fl k} = Jg^ih, k 2 ; Q) + jf^ih, k 2 ; Q). 



(4.18) 



For the evaluation of J'^ >1 (ki, k 2 \ Q) it is sufficient to calculate the trace of 
the Dirac matrices and integrate over k 



where J 2 (M N ) describes a divergent contribution depending on the cut-off 
Ad- Due to inequality Mn 3> Ad we can neglect J 2 (Mn) with respect to the 
convergent contribution. This corresponds too the accounting for the leading 
contributions in the large Nc expansion. Indeed, according to Eq.( |3.12D 
the contribution of divergent integrals is of order 0(1/ Nc) relative to the 
convergent ones. 

For the evaluation of J? 2 ^{k\,k 2 )Q) it is convenient, first, to integrate 
over k directions and then to calculate the trace over Dirac matrices. This 
gives 



J^{Kk2-,Q) = [l + 2J 2 (M^)][(Q + 2k 1 + k 2 ) a g 
+ {Q + 2k x + fc 2 ) V° + (Q + 2h + k 2 )»g af3 
-2(k 1 + k 2 ) a g^-2k^g'" a }, 



(4.19) 




xtr ( 7 /3 7 A 7 a 7^ + V + W 7 V) = 

= - ~[1 + 6 J 2 (M N )] [(3 Q + 2 h + fc 2 ) V + (3 Q + 2 h + k 2 fg 
9 

+ (3Q + 2h + k 2 ^g a % (- 



(4.20) 



Here we have used the integrals 



J n 2 i (M^-k 2 ) 4 
Summing up the contributions we obtain 




(4.21) 



J^(h,k 2 ;Q) = \ {Q a g^ + QPg» a + Q>*g a P) + % -[l + \j 2 {MK)} 

x [(2 h + k 2 ) a g^ + (2h + k 2 f g» a + (2 k x + k 2 y g a ?} 

+ [1 + 2 J 2 (M N )] [- 2 (k, + k 2 ) a gto - 2 k{ g^\. (4.22) 
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It is seen that the Q-dependence coincides with that obtained by means of 
the Gertsein-Jackiw method ( |4.17| ). Due to the arbitrariness of Q we can 
absorb by the Q-term the terms having the same Lorentz structure. This 
brings up the r.h.s. of ( |4.22j ) to the form 

+ [-2{k x + k 2 ) a g^ -2k^g^ a }. (4.23) 

Also we have dropped here the divergent contribution. This approximation 
is valid due to the inequality M N 3> A D and at leading order in the large N c 
expansion. 

The effective Lagrangian £Fi g .2c(a ; ) determined by the structure function 
( gggfl reads 



2 



£Fi g .2c(x) = ie^[$-a)d*Dl(x)D v (x)A v [ 



x 



6tt 2 lv ' ^ 

- (3 - a) Dfe) d u D u (x) A»(x) - bDfa) D v (x) A u (x) 

-(h-a) D\{x) D v {x) d u A^x) 

- (a - b) d u Dl(x) D»(x) A v {x) + b D^x) d u D"(a?) A v (x) 

+ 3 Dl(x) D u {x) {d" A v {x) - d v A»(x))]. (4.24) 

Due to the constraints d^D'Jx) = d^D^x) = some terms in the La- 
grangian ( |4.24j ) can be dropped out. This gives 

^Fig.2c( 2; ) = 

= ie^[-bDUx) DJx) d^A v (x) - (b - a) D]Xx) DJx) d v A"(x) 

- (a - b) d u Dl(x) D»(x) A v (x) + bD^x) d v D"{x) A v {x) 

+ 3 Dl(x) D u {x) {d"A v {x) - d v A»{x))}. (4.25) 

By using the relations d v D^ [x) = Dl (x)+d IJi Dl(x) and d u D^(x) = D v ^(x) + 
d ix D v {x) we can rewrite the Lagrangian ( fL25 ) as follows 

^Fig.2c(^) = 

= % e |^[- (a - b) Dt^x) A v {x) + bD^(x) A v {x) D^x) 

- (a - b) d v Dl(x) D u (x) A^x) + bD^x) d»D u (x) A v (x) 

- b D\{x) D v (x) d^A u {x) - (b - a) D\(x) D v (x) d v A»(x) 

+ 3 Dl(x) D u {x) {d»A v {x) - d v A»{x))}. (4.26) 
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The subsequent transformations we perform by applying the identity 

d u Dl(x)D u {x)A^\x) - Dl{x)d^D u {x)A u {x) = 

= Dl{x) D u {x) (d^A^x) - d v A»{x)) (4.27) 

being valid up to the contribution of a total divergence which can be omitted. 
Setting a = 2b we represent the effective Lagrangian (|4.27|) in the irreducible 
form 

£ F i g .2c0r) = ie^[bDl v {x)A»{x)D»{x)-bD^{x)A v {x)Dl{x)\ 

+ ie-^(-2b + 3)Dl(x)D v (x)F^(x), (4.28) 

where F^ u (x) = d^A u (x) — d u A tl {x) is the electromagnetic field strength 
tensor. Then, first two terms define the finite contributions to the renormal- 
ization constant of the wave function of the deuteron, whereas the last term 
coincides with the well-known phenomenological interaction Cc^{x) intro- 
duced by Corben and Schwinger [18] 

£ cs (x) = i e |^ (- 2 b + 3) D\{x) D v {x) F^(x) (4.29) 

for the description of the charged vector field coupled to an external electro- 
magnetic field. 

Thus, at leading order in the large Nc expansion the anomaly of the one- 
nucleon loop triangle VVV-diagram with vector (V) vertices defines fully the 
effective Lagrangian of the Corben-Schwinger type describing the deuteron 
coupled to an external electromagnetic field. In turn, the finite the contri- 
butions to the renormalization constant of the wave function of the deuteron 
we will identify below with those induced by a minimal inclusion of the elec- 
tromagnetic interaction: d^D v {x) — > (9 M + i eA v {x))D v (x). These terms 
are important for the correct definition of the effective Lagrangian of the 
deuteron coupled to an electromagnetic field. 

4.2 The phenomenological Aronson interaction 

The effective Lagrangian described by the diagram in Fig. 2d is defined 
by [15] 

J d 4 x C Fig . 2d (x) = J d 4 x J d S^^ 1 d frty 2 Dal3 ^ D t*( x i) A ^ x *) 
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a 2 1 

x e - ifcl ' xi e ~ ik2 - X2 e l{kl + k2 > x (- e) — j J a ^ vX {k u k 2 ; Q). (4.30) 



4tt 2 M 2 



In the structure function J a/3puX (ki,k 2 ;Q) is represented by the following 
momentum integral 



1 



d A k 

7T 2 i 



xtr < a 



a/3 



M^-k-Q 



■ a 



fl v 



1 



1 



'7 



M^-k-Q-ki M N -k-Q -k 1 -k 2 ) 

(4.31) 



The 4-vector Q = a k\ + b k 2 is an arbitrary shift of a virtual momentum, 
where a and 6 are arbitrary parameters. The Q-dependent part of the struc- 
ture function we obtain by using the Gertsein-Jackiw method [15] 



SJ ^\ kl ,k 2 ;Q) = 



1 

12 



tr( 7pC r a V^Q7 A + <r aP l P <?^l P Ql x + 



a/3 , 



1 



+ ( 7 a/3 Q<r^7 p7 V) = -tr((j Q/3 Q(7^7 A ). (4.32) 

6 



Now we should proceed to the evaluation of J al3 ^ vX {ki,k 2 ;Q). By analogy 
with J^ih, k 2 - Q) we get 



J a ^ x (h,k 2 ;Q) = 



r <¥k_ 

J 7T 2 i 



tr{a 



a/3 



M N + k + Q 
M?, - k 2 



1 + 



2k-Q 
M 2 -k 2 



a 



M N + k + Q + kt 
M 2 - k 2 



x 



x 



1 + 



1 + 



7 



2k- (Q + fci) 

M 2 - A; 2 
2 A; • (Q + h + k 2 ) 



x M N + k + Q + h + k 2 



M 2 -k 2 



M 2 -k 2 



d 4 k 



7T 2 i (M 2 - k 2 f 



+a a(3 a' x V(fc + Q + + fe)] + <7 a "(fc + QK"(fc + Q + £4)7 



tr {M^ [a a/3 (A; + QK V + a af3 a^(k + Q + hW 



a/3. 



xik + Q + h + k 2 )} 



1 + 



2 fc • (3Q + 2^ + k 2 ) 
M 2 -k 2 
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= / Pi (M* - P)3 tr {^nK^^ V + ° aP ^{Q + k)l X 
+a a/3 a^ 1 \Q + k + k 2 )\ - ^fcV^Q^V} 

+a a V^(3 Q + 2 fc x + £ 2 ) 7 A + a Q/ W(3 Q + 2 ^ + fc 2 )] 

-^A;V*' 3 (3Q + 2£; 1 + £ 2 KV} = 
o 

= J^ u \h, k 2 ; Q) + J^\h, k 2 ; Q). (4.33) 
Integrating over k we obtain 

J^ v \kiM,Q) = i[l + 2J 2 (M N )]tr(a^g^V) 

+ 1 tr [c^V^Q + ^) 7 A + <7 a V V(Q + fci + fe)], 

J^ uX (h,k 2 ;Q) = -^[-^ + J 2 (M N )]tr[a^(3Q + 2k 1 + k 2 )cr^ 1 x ] 

~tr [a a %3Q + 2k 1 + k 2 )v fMU l X 

+a a ?a^(3Q + 2k 1 + k 2 )l X + <7 Q/ VV(3Q + 2k x + k 2 )\. (4.34) 
Now we should sum up the contributions and collect like terms 

J a ^\h, k 2 - Q) = l - tr[^(Q - 1(2 h + k 2 ))a» V] 

+ I t r[ ( x a/3 ^^(fc 1 -fc 2 )7 A ] + -tr[<7 a/ VV(^i + 2jfe 2 )l. (4.35) 
6 6 

It is seen that the Q-dependence agrees with that obtained by the Gertsein- 

Jackiw method. Due to arbitrariness of Q the vector (2k\ + k 2 )/Q can be 

removed by the redefinition of Q. This gives 

J a ^ x (h, k 2 ; Q) = \ tr^Qa'V ) 
o 

+i tr [a^a^Ch - k 2 )l X ] + ^ tr [<r Q/ VV (^ + 2 fc 2 )]- (4.36) 

By evaluating the traces over Dirac matrices we obtain the structure function 
leading to the following effective Lagrangian [15] 

C Fig . 2d (x) = {-ie)^-^\lad x D^{x)D v »{x)A»{x) 



4tt 2 Ml 



3 
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+ -aD^»(x)d x D Xv A fl (x) + -(b + a)Dl v (x)D"\x)d^A x (x) 

+ 1 (6 - o) Dt„(x) D»\x) d x A"(x) - y Dl v {x) D» x (x) ^ A A (z) 

+ 8 D\ v {x) D v \x) (9" A x (x) - d x A"(x))] . (4.37) 

For the derivation of the effective Lagrangian ( 4.37|) we have used the equa- 
tion of motion 



d x D^{x) + d^D vX {x) + d v D^(x) = 0. 



(4.38) 



The analogous equation of motion is valid for the conjugated field. The term 
proportional to k\ contributing to the effective Lagrangian in the form of a 
divergence of the vector potential of the electromagnetic field d x A x (x) can be 
omitted by singling out the Lorentz gauge constraint for the electromagnetic 
potential, i.e. d x A x (x) = 0. 

Collecting like terms in ([4.37 ) we get 



-^Fi K .2d( a; ) 



le) 



4vr 2 M 2 



(b + a - 1) D^x) D vX {x) d» A x {x) 



+ ( 6 _ a - 3 ) DUx) D vX {x) d x A»(x) + -ad x D^ Xu (x) D vlx {x) ^{x) 



+ ~aD^(x)d x D x "{x)A tl {x) 



(4.39) 



The third and the fourth terms can be reduced by applying the equation of 
motion 

3 x D Xv (x) = -MlD v (x) 
and analogous for the conjugated field. Then, setting 



b + a - 1 



-b + a + 3 



(4.40) 



we obtain 6 = 2 and bring up the effective Lagrangian ( |4.39|) to the following 
irreducible form 



^Fig.2d(^) 

2g* 



i e 



+ ie 



3tt 2 
2g\ 1 



a [- Dl„{x) A v {x) D»{x) + D^{x) A v (x) D^x)} 



3tt 2 M 2 



(1 + a)Dl v {x)D vX (x){d x A»(x)-d»A x {x)). (4.41) 
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The first two terms define the finite contributions to the renormalization 
constant of the wave function of the deuteron, whereas the last term coin- 
cides with the well-known phenomeno logical interaction Ca{x) introduced 
by Aronson [19] 

C A (x) = ie M _L (i + a ) Dl {x ) D uX (x) F x »(x) (4.42) 

for the description of the charged vector field coupled to an electromagnetic 
field. 

Thus, we have shown that the anomaly of the one-nucleon loop triangle 
V^TT-diagram, where V and T stand for the vector and tensor vertices deter- 
mined by the Dirac matrices 7° and a^, respectively, calculated at leading 
order in the large Nc expansion defines fully the phenomenological Aronson 
Lagrangian describing the deuteron coupled to an external electromagnetic 
field. 

4.3 The magnetic dipole and electric quadrupole mo- 
ments of the deuteron 

The effective Lagrangian describing both the magnetic dipole and electric 
quadrupole moments of the deuteron is determined by the sum of £Fi g .2c(^) 
and ^Fig.2d(^) given by Eq. (|4.28|) and Eq.( |4.41| ), respectively, and reads 



t (rA A v !nA DM/ 



^^_4o^^ 



+ ie-^{-2b+ 3)DUx)D u (x)F^{x) 

+ ie (l + a)^^Dl(x)D»\x)F x »(x), (4.43) 

where a and b are arbitrary parameters related to the ambiguities of the one- 
nucleon loop diagrams with respect to a shift of a virtual nucleon momentum. 
We consider them as free parameters of the approach [15]. 

In order to fix these parameters it is convenient to write down the total 
effective Lagrangian of the physical deuteron coupled to an external electro- 
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magnetic field 



x 



+ ie 



~Dl(x)D^(x)+M^Dl(x)D' , 
bg2v ;* ag2T Dl(x)A»(x)D»(x) 



x 



i e ■ 



6tt 2 



D 



a; 



+ ie 



+ 



6vr 2 v 
i e (1 + a 



26 + 
' 3tt 2 



3) D„(z) 
1 



M 2 



t/A/ 



(4.44) 



Two terms having the structure D^ v (x) A^(x) Dl(x) and Dj lu (x)A' 1 (x)D u (x) 
should describe the interaction of the deuteron with an external electromag- 
netic field included by a minimal way, whilst the last two terms are respon- 
sible for the non-trivial contributions to the magnetic dipole and electric 
quadrupole moments of the deuteron. In terms of the parameters of the ef- 



fective interactions Eq. fl4.44p the magnetic dipole moment /id, measured in 
nuclear magnetons, and the electric quadrupole moment Qd, measure in fm 2 , 
of the deuteron are given by 



Qd 



[l + a 



9t 



3tt 2 



[2 + 2a) 



+ (3-26) 
(3 



9t 



3tt 



127T 2 ' 

2b) & 
6tc 2 



Ml 



at the constraint 



6vr 2 



- 2a 



9t 
3vr 2 



1 



(4.45) 



(4.46) 



reducing the first two terms in effective Lagrangian Eq. (|4.43 ) to the standard 
minimal form which can be obtained from the effective Lagrangian of the free 
deuteron field by the shift d^D u (x) — > (d^ + i e A^{x))D v {x). 

Retaining the former relation between the electric quadrupole moment 
and the coupling constant gy, Qd = 2g Y /7i 2 Ml [15] that gives gy = 11.319, 
we express the parameters a and b in terms of the coupling constants gy, g^ 
and the magnetic dipole moment /io: 

- 3 7T 2 / *f 
- — + ^ 



2g$ 



9 

b=- 
2 



37T 2 

9y 



/i D . 



(4.47) 
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Substituting Eq.( [4.47| ) in Eq.( |4.46| ) we get the relation between coupling 
constants gy, and the magnetic dipole moment /id 



9t 



^^ + ^7T 2 (l + ^ D ] = 0.7995V. (4.48) 




The numerical value g T = 0.799 gy we obtain at gy = 11.319 [15], /xd = 0.857 
[22] and Nc = 3. The sign of the coupling constant #t should coincide with 
the sign of the coupling constant gy. For the opposite sign the binding energy 
£d given by Eq. (|4.8|) becomes negative that means the absence of the bound 
neutron-proton state with quantum numbers of the deuteron. 

For the evaluation of the binding energy e~d determined by Eq. ( |4.8| ) we 
should keep only the leading contribution to the coupling constant g? in the 
large Nc expansion, i.e. 




9r = \ o9v + 0(l/JN c ). (4.49) 



Substituting this relation into Eq.( |4.8|) we can describe the experimental 
value of the binding energy of the deuteron Ed = 2.225 MeV at the cut-off 
Ad = 46.172 MeV. The spatial region of virtual nucleon field fluctuations 
forming the physical deuteron related to this value of the cut-off 1/Ad ~ 
ro = 4.274 fm agrees good with the experimental value of the radius of the 
deuteron = (4.31895 ± 0.00009) fm [22]. This result confirms estimates 
obtained in Ref. [15]. 

The effective Lagrangian of the deuteron field coupled to an external 
electromagnetic field is given by 

= ~\ [(«/. - ieA^DUx) - [d v - ieA v {x))D\{x)\ 
x [(&* + % e A"(x))D u (x) - (d u + i e A v {x))D*{x)\ + M£ (x)D"(x) 

+ ieL D - ^QnM^J D^x) D v {x) F^{x) 



+ ie U D + ~Q»M*\ w Dl(x)D»\x)F x »(x). (4.50) 



1 „ „o\ 1_ 

The term of order 0(e 2 ) can be also derived in the NNJL model by using 
shift ambiguities of one-nucleon loop diagrams. This term is required by the 
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electromagnetic gauge invariance of the effective Lagrangian of the deuteron 
field coupled to an external electromagnetic field, but it does not affect on 
the electromagnetic parameters of the deuteron which are of order 0(e). 

5 Conclusion 

We have shown that the Nambu-Jona-Lasinio model of light nuclei or 
the NNJL model as well as the ENJL model with chiral £7(3) x £7(3) sym- 
metry [4-11] is motivated by QCD. The NNJL model describes low-energy 
nuclear forces in the nuclear phase of QCD in terms of one-nucleon loop ex- 
changes. One-nucleon loop exchanges provide a minimal way of the transfer 
of nucleon flavours from an initial to a final nuclear state and allow to take 
into account contributions of nucleon-loop anomalies. These anomalies are 
related to high-energy fluctuations of virtual nucleon fields, i.e. the NN 
fluctuations, and fully determined by one-nucleon loop diagrams [29-31]. 
The dominance of contributions of one-nucleon loop anomalies to effective 
Lagrangians describing low-energy interactions of the deuteron coupled to 
itself, nucleons and other particles we justify within the large Nq< expansion 
in QCD with SU(Nc) gauge group at Nq — > oo. It is well-known that 
anomalies of quark-loop diagrams play an important role for the correct de- 
scription of strong low-energy interactions of low-lying hadrons [4-11]. We 
argue an important role of nucleon-loop anomalies for the correct description 
of low-energy nuclear forces in the nuclear physics. 

It should be emphasized that nucleon-loop anomalies can be interpreted 
as non-trivial contributions of the non-perturbative quantum vacuum - the 
nucleon Dirac sea [32]. In nuclear physics the influence of the nucleon Dirac 
sea on low-energy properties of finite nuclei has been analysed within quan- 
tum field theoretic approaches in the one-nucleon loop approximation [33]. 
Unfortunately, in these approaches contributions of one-nucleon loop anoma- 
lies have not been taken into account. The NNJL model allows to fill this 
blank. 

For the derivation of the NNJL model from the first principles of QCD 
we distinguish three non-perturbative phases of QCD: 1) the low-energy 
quark-gluon phase (low-energy QCD), 2) the hadronic phase and 3) the nu- 
clear phase. Skipping over the intermediate low-energy quark-gluon phase by 
means of the integration over high- and low-energy quark and gluon fluctua- 
tions one arrives at the hadronic phase of QCD containing only local hadron 
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fields with quantum numbers of mesons and baryons coupled at energies be- 
low the SB%S scale A x ~ 1 GeV. The couplings of low-lying mesons with 
masses less than the SB%S scale to low-lying octet and decuplet of baryons 
can be described by Effective Chiral Lagrangians with chiral U(3) x U(3) 
symmetry. 

Integrating in the hadronic phase of QCD over heavy hadron degrees of 
freedom with masses exceeding the SB%S scale one arrives at the nuclear 
phase of QCD which characterizes itself by the appearance of bound nucleon 
states - nuclei. At low energies the result of integration over heavy hadron 
degrees of freedom can be represented in the form of phenomenological local 
many-nucleon interactions. Some of these interactions are responsible for 
creation of many-nucleon collective excitations which acquire the properties 
of observed nuclei through nucleon-loop and low-lying meson exchanges. 
This effective field theory describes nuclei and processes of their low-energy 
interactions by considering nuclei as elementary particles represented by local 
interpolating fields. 

Following this scenario of the description of nuclei and their low-energy 
interactions from the first principles of QCD the deuteron should be produced 
in the nuclear phase of QCD by a phenomenological local four-nucleon in- 
teraction as the Cooper np-pair with quantum numbers of the deuteron. 
The low-energy parameters of the physical deuteron, i.e. the binding en- 
ergy, the magnetic dipole /id and electric quadrupole Qd moments and so, 
the Cooper np-pair acquires through one-nucleon loop exchanges. We have 
shown that the main part of the kinetic term of the effective Lagrangian of 
the free physical deuteron field is induced by the contribution of high-energy 
(short-distance) fluctuations of virtual nucleon fields related to the anomaly 
of the one-nucleon loop VV-diagram with two vector vertices. 

In turn, the magnetic dipole ^ and electric quadrupole Qd moments of 
the physical deuteron are fully determined by high-energy (short-distance) 
fluctuations of virtual nucleon fields related to the anomalies of the triangle 
one-nucleon loop VVV- and VTT-diagrams. Thus, high-energy (short- 
distance) fluctuations of virtual nucleon fields related to anomalies of one- 
nucleon loop diagrams play a dominant role for the correct description of 
electromagnetic properties of the physical deuteron in the NNJL model. 

As regards low-energy (long-distance) fluctuations of virtual nucleon 
fields they give a significant contribution only to the binding energy of the 
deuteron e-Q. The strength of low-energy (long-distance) fluctuations of vir- 
tual nucleon fields is restricted by the cut-off A D = 46.172 MeV. The spatial 
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region of virtual nucleoli field fluctuations forming the physical deuteron re- 
lated to this value of the cut-off 1/A D ~ r D = 4.274 fm agrees good with the 
experimental value of the radius of the deuteron td = (4.31895 ± 0.00009) fm 
[22]. This confirms our estimates obtained in Ref. [15]. 

It is well-known that in the potential model approach to the description of 
the deuteron the electric quadrupole moment of the deuteron Q D is caused by 
nuclear tensor forces which are of great deal of importance for the existence 
of the deuteron bound np-state [34]. 

The proportionality of the coupling constant of the phenomenological 
local four-nucleon interaction Eq.( |3.3| ), responsible for creation of the Cooper 
np-pair with quantum numbers of the deuteron, and the binding energy of 
the deuteron e D Eq. Q3.20D to the electric quadrupole moment Qd testifies an 
important role of nuclear tensor forces for the formation of the deuteron in 
the NNJL model. 

To the evaluation of one-nucleon loop diagrams defining effective La- 
grangians describing processes of low-energy interactions of the deuteron 
coupled to itself and an electromagnetic field we apply expansions in powers 
of the momenta of interacting particles and keep only leading terms of the 
expansions. This approximation can be justified in the large Nc expansion. 
Indeed, in QCD with the SU (Nc) gauge group at Nc — > oo the nucleon mass 
is proportional to the number of quark colours [17]: Mn ~ Nc- Since for the 
derivation of effective Lagrangians describing the deuteron and amplitudes 
of low-energy nuclear processes all external momenta of interacting parti- 
cles should be kept off-mass shell, the masses of virtual nucleon fields are 
larger compared with the external momenta. An expansion of one-nucleon 
loop diagrams in powers of 1/Mn giving an external momentum expansion 
corresponds to the expansion in powers of 1/Nq. In this case the leading 
order in the large Nc expansion gives the leading order contributions in the 
expansion in powers of external momenta of interacting particles. We should 
emphasize that anomalous contributions of one-nucleon loop diagrams are 
determined by the least powers in external momentum expansions. Thereby, 
the dominance of contributions of nucleon-loop anomalies to effective La- 
grangians describing low-energy nuclear forces in the NNJL model is fully 
supported by the large Nc expansion. The accuracy of this approximation is 
rather high. Indeed, the real parameter of the expansion of one-nucleon loop 
diagrams is 1/M^ ~ 1/N C but not 1/Mn ~ 1/Nc- Thereby, next-to-leading 
corrections should be of order 0(1/N C ). 

The inclusion of the interaction of the deuteron field with the tensor 
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nucleoli current Eq. ([4.1|) has given a possibility of the self-consistent descrip- 
tion of the electromagnetic properties of the deuteron, the magnetic dipole 
moment /xd and the electric quadrupole moment Qn, in terms of effective 
interactions of the Corben-Schwinger and Aronson kinds induced by one- 
nucleon loop diagrams. By fitting the experimental values of the magnetic 
dipole moment /i D = 0.857, measured in nucleon magnetons /in = e/2M N , 
and the electric quadrupole moment Q D = 0.286, measured in fm 2 , supple- 
mented by the requirement of the electromagnetic gauge invariance of the 
effective Lagrangian of the deuteron field coupled to an external electromag- 
netic field we have got the relation between the coupling constants gy and 
9t- 9t — 0.799 gy calculated at N c = 3. At leading order in the large N c 
expansion we get g T = ^3/8 gy + 0(l/y/N c ). This relation agrees good 
with that obtained in Ref. [14] (see Eq.(16) of Ref. [14]). Due to this rela- 
tion the experimental value of the binding energy of the deuteron can be 
described by the cut-off A D = 46.172 MeV. This corresponds to the spa- 
tial region of virtual nucleon field fluctuations forming the physical deuteron 
1/Ad ~ td = 4.274 fm agreeing good with the experimental value of the 
radius of the deuteron r D = (4.31895 ± 0.00009) fm [22]. 

For further applications of the NNJL model to the description of low- 
energy nuclear reactions of astrophysical interest we anticipate the results in 
agreement with those obtained in Refs. [35,36]. 

The quantum field theoretic scenario to treating nuclei as many-nucleon 
collective excitations induced by phenomenological local many-nucleon in- 
teractions allows a plain extension of the NNJL model by the inclusion of 
light nuclei 3 He, 3 H and 4 He as three- and four-nucleon collective excitations. 
The binding energies and other low-energy parameters of these excitations 
should be determined through nucleon-loop and low-lying meson exchanges. 

On this way it is important to notice that the spinorial structure of 
the operators of three-nucleon densities coupled to the 3 He and the 3 H is 
very much restricted. One can show that only the three-nucleon densities 
\p c (x) r y tJ,n f 5 p(x)] r y tJl ,n(x) and [n c (x)'y^'y 5 n(x)} r y fl p(x) can lead to the appearance 
of the bound 3 He and 3 H state, respectively. At the quantum field theoretic 
level this result explains a well-known experimental fact of the compensation 
of spins and magnetic dipole moments of pp and nn pairs inside nuclei which 
has been put into the foundation of the shell-model of nuclei [34]. 

The extension of the NNJL model by the inclusion of 3 He, 3 H and 4 He 
would give a possibility to analyse within the NNJL model the reactions of 
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the p-p chain [37] started with the reaction p + p^D + e + + z/ c and 
to apply the extended version of the NNJL model to the description of the 
reactions p + D — > 3 He + 7, p + 3 He — > 4 He + e + + v e and so on. 

Chiral perturbation theory can be naturally incorporated into the NNJL 
model [35] in terms of Effective Chiral Lagrangians with chiral (7(3) x (7 (3) 
symmetry describing low-lying baryons and mesons interacting at low ener- 
gies [4-11]. 

The quantum field theoretic description of the deuteron within the NNJL 
model can be also of use for the analysis of the properties of dibaryons. 
Indeed, following Oakes [38] the deuteron can be considered as a component 
of the 5'(7(3) flavour decuplet 10/ of dibaryons with Y = 2 and 1 = 0, where 
Y and I are the hypercharge and the isotopical spin, respectively. In the 
chiral limit the binding energies, the magnetic dipole and electric quadrupole 
moments of dibaryons of the decuplet 10/ should be equal. The splitting of 
the parameters of the components of the decuplet 10/ can be obtained within 
Chiral perturbation theory incorporated into the NNJL model. 
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Figure caption 



• Fig.l. One-nucleon loop diagrams contributing in the NNJL model to 
the binding energy of the physical deuteron, where rf = C n T is the 
field of anti-nucleon. 

• Fig. 2. One-nucleon loop diagrams describing in the NNJL model the 
effective Lagrangian of the deuteron coupled to an electromagnetic field 
through the magnetic dipole and electric quadrupole moments, where 
n c = C n T is the field of anti-nucleon. 
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